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1. Introduction 


In the last few years, several new discoveries have shed light on the AdS/CFT correspon¬ 
dence i I I- This correspondence maps strings moving in an AdS^ x background 
to an AA = 4 supersymmetric Yang-Mills (SYM) theory. The eigenvalues of the dilatation 
operator are mapped to the energies of closed string states Q. A step in understanding 
this duality better was the discovery that the dilatation operator of the AA = 4 SYM is 
proportional to the Hamiltonian of an integrable spin chain |,i0- 

Recently, progress has been made to extend the gauge/gravity-correspondence, in con¬ 
text of spin chains, towards more realistic models with less supersymmetry 
For instance, if the background geometry for the string is AdS^ x W, where W is some 
compact manifold, the dual gauge theory should still be conformal. Other geometries, 
mainly orbifolds of AdS^ x corresponding to non-conformal theories have also been 
investigated | 
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The success in using spin chains to study the duality beyond the BMN limit motivates 
studies of integrability of deformed correlators. One question that naturally arises in this 
context is whether integrability is related to supersymmetry, conformal invariance or have 
more geometrical reasons. 
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The Leigh-Strassler deformations [^] preserve AA = 1 supersymmetry and conformal 
invariance, at least up to one loop. It is hence of great interest to investigate if there exist 
points in the parameter space where the dilatation operator is mapped to an integrable 
spin-chain Hamiltonian. This question has been under investigation in |^, . In 

|n] this deformation was studied in a special case corresponding to a g-deformed (often 
called /3-deformed) commutator. It was found that for the sector with three chiral fields 
the dilatation operator is integrable for q equals root of unity. 

In reference p0| , a way of generating supergravity duals to the /3-deformed field theory 
was introduced, and in &I 0 agreement between the supergravity sigma model and 
the coherent state action coming from the spin chain describing the /3-deformed dilatation 
operator was demonstrated. This way of creating supergravity duals was used in |Q to 
construct a three-parameter generalization of the /3-deformed theory. The gauge theory 
dual to this supergravity solution was found in [^, |T^ for q = with real, corre¬ 
sponding to certain phase deformations in the Lagrangian. This gauge theory is referred 
to as twisted SYM, from which the /3-deformed theory is obtained when all the 'jj = (3. 
The result is that the theory is integrable for any q = with 7 j real |18]. The general 
case with complex 7 ^ is not integrable MM- 

In the present work, the g'-deformed analysis is extended to the more general Leigh- 
Strassler deformations with an extra complex parameter h, in order to find new integrable 
theories. A site dependent transformation is found which relates the 7 j-deformed case to a 
site dependent spin-chain Hamiltonian with nearest-neighbour interactions. In particular 
when all are equal, the transformation relates the g-deformed theory to the /i-deformed 
theory, i.e. the theory only involving the parameter h. In particular, we find a new 
R-matrix, at least in the context of AA = 4 SYM, for g = 0 and h = e*® with 6 real. 
We also find all R-matrices with a linear dependence on the spectral parameter which 
give the dilatation operator. A general ansatz for the R-matrix is given. Unfortunately, 
the most general solution is not found. However, we find a new hyperbolic R-matrix 
which corresponds to a basis-transformed Hamiltonian with only diagonal entries |^|. A 
reformulation of the general R-matrix shows that the structure of the equations obtained 
from the Yang-Baxter equations resemble the equations obtained in the eight vertex model. 
This gives a clear hint how to proceed. 

In the dual supergravity theory, some attempts to construct backgrounds for non-zero 
h have been done I 


. Apart from the five-flux there is also a three-flux. A step going 
beyond supergravity was taken in where the BMN limit was considered. We hope 
our results will make it easier to find the supergravity dual of the general Leigh-Strassler 
deformed theory. 


2. Marginal deformations of N'=4 snpersymmetric Yang-Mills 

To study marginal deformations of AA = 4 SYM with SU{N) gauge group, it is convenient 
to use M = 1 SYM superfields. The six real scalar fields of the AA = 4 vector multiplet 
are combined into the lowest order terms of three complex M = 1 chiral superfields 4>o, 4>i 
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and d> 2 . It is well known that the N = 1 superpotential 


wm=i = , 


( 2 . 1 ) 


where is the coupling constant, describes a finite theory at one-loop if the following 

two conditions are fulfilled ( 2 ^, 

3C2{G) =Y,T{Ai), and 


Ci^dCf^L = 2g^TiAi)5j>5^j 


( 2 . 2 ) 


The constant 6*2 (G) is the quadratic Casimir operator defined here^ as 6*2 (G) • 1 = 6ahT%T\ 
where A is the adjoint representation of the group G which in the present context is the 
symmetry group SU{N). The constant T{M) is defined through T{M)6°'^ = Tr(T^T^) for 
the representation M. The first condition of (p.2|) implies that the /3-function is zero. For 


an SU{N) group with the superpotential (^) this is automatically fulfilled. The choice 
^abc^ = fabc therefore gives a super conformal Af = I theory at one-loop. However, 

there are more general superpotentials satisfying the one-loop finiteness conditions. To ex¬ 
plore marginal deformations of AA = 4 SYM we consider the Leigh-Strassler superpotential 

IT! 


W = ^Ae^^^Tr [[$/, $j] ^k] + [{$,, $j} <^k] 


(2.3) 


where is totaly symmetric. The coupling constants can now be written as = 

^^IJK fabc _|_ ■^^Tr [{Ta, Tft} Tc]. The non-zero couplings are chosen to be = A(1 — 

q)/{l + q) and = 2Xh/{l + q)- In terms of the deformation parameters q and h the 
superpotential (|2.3|) becomes 


W = 


2A 

1 + g 


Tr 


$o4>i$ 2 - ^ ($0 -h -h ^ 2 ) 

O 


(2.4) 


This deformed superpotential will be our main focus. 

The presence of q and h in the superpotential (^~i ) breaks the SU (3) symmetry in the 
chiral sector. What is left of the symmetry is a Z 3 x Z 3 symmetry. The hrst Z^, permutes 
the ‘h’s and the second takes ‘ho —> w4>o, 4)i ^ a;^<I>i and ^2 —> ‘^ 2 , where a; is a third root 
of unity. 

The one-loop finiteness condition (^) is satisfied if 


g 


2 


A2 

(l + g )2 




(2.5) 


In the large-iV limit, which we consider, the relation ( p.5| ) becomes even more simple. The 
one-loop finiteness condition ( 2 . 2 ) also implies that the scalar field self-energy contribution 
from the fermion loop is the same as in the AA = 4 scenario, due to the fact that the fermion 
loop has the contraction The parameters in (^) span a space within which 


^Our conventions are: r“ are the SU{N) group generators, satisfying T°‘ = r“h The normalization of 
r“ is given by Tr(T“r'’) = 5“*’/2 from where it follows that Tt:{T%T\) = N5°'*’ in the adjoint representation 

T. 
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there exists a manifold, or perhaps just isolated points, (5{g, A, q,h) = 0 of superconformal 
theories to all loops |15]. In the limit g —> 1 and /i —> 0 the AA = 4 SYM is restored. 
Marginal deformations away from this fixed point will be explored in the following sections 
by means of integrable spin chains. 


3. Dilatation operator 


From the Leigh-Strassler deformation (^) of the AA = 4 SYM theory it is possible to 
obtain the dilatation operator in the chiral sector. In this sector, the only contribution 
is coming from the F-term in the Lagrangian, under the assumption that the one-loop 
finiteness condition (|2.2| ) is fulfilled. The scalar field part of the F-term can be expressed 
in terms of the superpotential as 


Cf 


<9IF 

2 

dW 

2 

dW 

9$o 

-b 

d<^i 

+ 

0^2 


(3.1) 


Using (pQ, (pi and (p 2 to denote the complex component fields, the Lagrangian becomes 
(omitting the overall factor 2\/[1 + q) in (2.4)) 


Cf = Tr [cpicpi+icpi+icpi - qcpi+icpicpi+icpi - q*cpicpi+icpicpi+i] 

+ Tr [qq*(pi+icpi^i^i+i - qh*cpi+icpi^i+ 2 ^i +2 - Q*h(pi+ 2 (pi+ 2 (pi(pi+i] 

T Tr ^h<piF24^i+2p^i+lP^i T h (pi(piJf-l(pi-\-2(pi+2 T hh (pifpifpifpi^ , (3.2) 


where a summation over i = 0,1,2 is implicitly understood and the indices of the fields (pi 
are identified modulo three. In order to see how the dilatation operator acts on a general 
operator O = ■■■(piL loop order in the planar limit we calculate the 

Feynman graphs and regularize in accordance with [IT, Dj. The vector space, spanning 
these operators, can be mapped to the vector space of a spin-1 chain (see m for details). 
We define the basis states |0), |1) and |2) for the spin chain which correspond to the 
fields (po^ (pi and (p 2 - By introducing the operators Eij, which act on the basis states as 
Eij\k) = 5jk\i), the dilatation operator can be written as a spin-chain Hamiltonian with 
nearest-neighbour interactions, i.e. A = where 




-qEi,,,E’p^l,-(fEl,,iE^ 


/ + 1 


/ + ! 


'/+! u jpl 




+ + h*EE^2E'[tl,+2 + hh*EEEl 


./ + 1 


l.i 


(3.3) 


The direct product between the operators E^j is suppressed. If we use the convention 



/l\ 


(“) 

/o\ 

|0) = 

0 

|i) = 

1 |2) = 

0 


[v 


[«/ 



(3.4) 
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the Hamiltonian can be expressed as the matrix 


/ h*h 

0 

0 

0 

0 

h 

0 

* 

1 

0 \ 

0 

1 

0 

- q * 

0 

0 

0 

0 

h* 

0 

0 

q*q 

0 

-qh* 

0 

-q 

0 

0 

0 

-q 

0 

q*q 

0 

0 

0 

0 

-qh* 

0 

0 

1 

* 

0 

h*h 

0 

h 

0 

0 

h* 

0 

0 

0 

0 

1 

0 

* 

-q 

0 

0 

0 

-q 

0 

h* 

0 

1 

0 

0 

1 

* 

0 

0 

0 

0 

-q 

0 

q q 

0 

V 0 

h 

0 

1 

* 

0 

0 

0 

0 

h*h / 


(3.5) 


|. The existence of a 
a), is crucial for the 


We will now search for special values of the parameters h and q for which the spin- 
chain Hamiltonian (3.3) is integrable. When h is absent, the analysis simplifies consider¬ 
ably, because the usual S-matrix techniques can be used 00 
homogeneous eigenstate, an eigenstate of the form | a)( 8 ) | a)... <5 
S-matrix techniques to work. From this reference state, excitations can be defined. In this 
context, the state |a) is a pure state, that is, one of the states | 0 ), | 1 ) or | 2 ). 

When h is non-zero, the analysis become significantly harder. The only values for the 
parameters, for which it is possible to define a homogeneous eigenstate are q = 1 + e®27rn/3^ 
or g = —1 and h = where n is an arbitrary integer. In these cases the homogeneous 

eigenstates are 


, . i27rm 

=| 0 ) + e— 


1 ) + e“ 


I 2), 


m G Z. 


(3.6) 


Clearly, the two Z 3 symmetries are manifest. For q = 1 + the eigenvalues are 

zero, thus the corresponding states are protected. This case is related to the q-deformed 
Hamiltonian by a simple change of variables. We introduce a new basis 


| 0 ) 
| 1 ) 
12 ) 


i27m 

e 3 


(|0) + |i) + 2)). 


^(| 0 ) + e-|l) + e--| 2 )), 

1 , i2'K i27r . 

^(| 0 ) + e--|l) + e-| 2 )), 


(3.7) 


where n is an integer. It will shortly be shown that the phase shift in |0) will imply that 
a phase can be transformed away from h. The Hamiltonian expressed in the new 

basis (^) takes the same form as (^), but with new parameters q and h and an overall 
proportionality factor 

't27r _ i27r _ i27rrt 

e 3 — qe 3 + he 3 . (3.8) 

The new parameters q and h can then be expressed in terms of the old parameters as 


q = 

h = 


i27r —i27r _ 227rn 

ge 3 — e~^ — he~ 3 


i27T —i27r _ i27rn 

e 3 — qe 3 + he 3 


_ i2'Kn 

1 — q + he 3 

^27r —i2'n _ i27rn 

e 3 — ge 3 + he 3 


(3.9) 

(3.10) 
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The case q = + 1 corresponds to the q-deformed case, and if /I’s imaginary 

part comes from the phase e®27rn/3^ remaining part is phase independent. This is in 
agreement with reference |29|. The integrable case q equals a phase will correspond to the 
case q = /ie-*27rn/3 _j_ h = pe*27rn/3 p and q being real. It is also clear that the 

case q = —1 and h = is related by the change of basis to a Hamiltonian of the form 


u \ ^ Q r F* + l _L F* + l _L F* r?' 

M — 2 _^ ^ [£/ 22n/22 + -^ 00-^00 + -^ 11-^11 


'i+ll 


(3.11) 


This case looks perhaps trivial, but it is not. The different eigenvalues equal 3n with 
n = 0,1, 2,... , L — 2, L. Note that the value L — 1 is absent for this periodic spin chain^. 
The states have a large degree of degeneration. 

For other values of q and h, a reference state does not have a precise meaning. Hence, 
we cannot adapt the S-matrix formalism. Instead, we will try to find an R-matrix, from 
which the Hamiltonian (3.3) is obtainable. The existence of an R-matrix R{u), depending 
on the spectral parameter u, is sufficient to ensure integrability. All R-matrices necessarily 
have to satisfy the Yang-Baxter equation 

Ru{u - v)Ri3{u)R23{v) = R23iv)Ri3{u)Ri2iu - v). (3.12) 

The Hamiltonian can be obtained from the R-matrix through the following relation 

(3.13) 


d 


P—R(u)\u=uo = H, 

where V is the permutation operator, with the additional requirement R{uq) = V for some 
point u = uq. 


4. A first look for integrability 


In this section, we will show how the transformation of basis (3.7) combined with a position 
dependent phase shift, sometimes called a twist, gives rise to new interesting cases of 
integrability. In |]^], the (7-deformed case was studied. It was shown that for q equals a root 
of unity, the phases can be transformed away into the boundary conditions. Furthermore, 
it was shown in |]^ that the integrability properties do not get affected for any q = 
where f3 is real. It was also established that a generalised SYM Lagrangian deformed with 
three phases 7^, instead of just one variable, is integrable. The deformed theory is referred 
to as the twisted (or 7-deformed) SYM and the corresponding one-loop dilatation operator 
in the three scalar sector is 

_ jpl I jpl l?^+l I jpl l?^+l 

H - -t/Qo-t/ii + + -^22-^00 

p*71 fpl Tpl Tpl+^ p*73 fp/ 

e -f- e 1 :^ 21^12 ^ ^ -^ 02-^20 

I 


+ e-^^^E[2Eif + e-^^^E^2oEo2 


iZ+l 


+ 


Tpl Tpl-\-^ I Tpl Tpl-\-^ I Tpl IT’^+l 


-^00-^22 


(4.1) 


^Excitations can be created if two states of the same number are next to each other. For example the 
state |112012) has energy three and the next highest energy state is |111112) with energy 4x3. The state 
with the highest energy, equals to 6 x 3, is 
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A natural question to ask is if the phases can also be transformed away in a generic 


Hamiltonian of the form (3.3). If both q and h are present we can not, at least in any simple 
way, transform away the phase of the complex variables. However, when q = it is 

possible to do a position dependent coordinate transformation 


|6)fc = e*2-/3|o),, |i)fc = e'2W3|i)^^ and |2)fc = 


(4.2) 


as in so that the phase of q is transformed away. Here, k refers to the site of the 
spin-chain state. This transformation changes the generators in the Hamiltonian as 


F' = c" 


-e: 


n,n+m * 


(4.3) 


This kind of transformation of basis generally results in twisted boundary conditions. Thus, 
the periodic boundary condition |o)o = |a )l for the original basis becomes in the new basis 


|6)o = |6 )l 


. " 1 l^TTL, , ~ , 

|l)o = e—|1 )l. 


and 


, ~ , —i 27 rL , ~ , 

|2)o = e—|2)i, 


(4.4) 


where L is the length of the spin chain. A consequence is that the system is invariant under a 
rotation of q by introducing appropriate twisted boundary conditions (4.4). As an example, 
the q-deformed Hamiltonian with periodic boundary conditions with q = + 1 

(see text above ( 3.11 )), is equivalent to -|- 1 with twisted boundary 

conditions. Hence, the following cases are integrable 


h = pe 


X^TT/l 

q={l+p)e 3 


and 


i27rm 

q = -e 3 , 


l27T7l 

h = e 3 


(4.5) 


where p is real and can take both negative and positive values and n and m are arbitrary 
independent integers. 

One can actually combine the twist transformation above with the shift of basis (13) 
in a non-trivial way. This combination will turn out to give a relation which maps the 
Hamiltonian with arbitrary q and vanishing h into the Hamiltonian with vanishing q and 
arbitrary h. The periodic boundary condition will, however, change for spin chains where 
the length is not a multiple of three. 

In terms of matrices the transformation can be represented as follows. Let us represent 
the shift of basis (3.7) by the matrix T (with n set to zero) 


T = 


1 


v/3 


11 1 \ 

I g*27r/3 ^-i2TTl3 

^ ^-i2irl3 gi27r/3 


and the transformation matrix related to the phase shift (| 
shift in the zero state |0)) 

/l 0 0 


(4.6) 


by (but without the phase- 


Uk = 


0 0 
0 0 


(4.7) 


®Note that the phase factor in |0) is not position dependent, it was only added in order to cancel the 
extra phase which would have appeared in front of the terms having h in them. 
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The transformation that takes the g-deformed to the /i-deformed Hamiltonian is then 


where 


H = TiHTf 1, 


Ti = (T ® T){Uk ® ® T-i) 


(4.8) 


(4.9) 


Acting with this transformation on the Hamiltonian (|3.3|) we get the new Hamiltonian 


jpl ZT'^+l I 


:a ’ 


Up to an overall factor, the transformation (^) change the couplings as 

q ^ 0 and h = 0 g = 0 and h = —l/q 

In terms of states, the map ([4.8|) generates the following change 


(4.10) 


(4.11) 


l“)l+3fc ^ l« — l)l+3fe ) |“)2+3fc ^ l« + l)2+3fe ) |“)3fe ^ l“)3fe ) 


(4.12) 


where a takes the values 0, 1 or 2. Let us investigate how the transformation ( |4.10| ) affect 
the boundary conditions. Prom equation ( 4.12|) we see that the original periodic boundary 
conditions |a)o = \a)L translate into 


| 0 " 


\^new)^^\^new)^ and 


= lU 


(4.13) 


if the length L of the spin chain is one modulo three and the opposite, |0"'®'^)o = etc, 

for the two modulo three case. If the length is a multiple of three the boundary conditions 
remain the same. 

If we start from the Hamiltonian of the y-deformed SYM (4.1), the transformation 
r^) leads to the Hamiltonian 


W 


i,i+i _ 


jpl Tpl+l I Tpi Tpl+^ I Tpi ir^+l 
-^ 00-^11 + -^ 11-^22 + ^22^00 

-i'ys-l pi IT'^+I I p-i'Ti-i pi pl+^ I p-il 2 -l pi K ^+1 

e -C'20-^21 + ® -^01-^02 + ^ -^12-^10 

A'i^-‘eLeA^ 


_L £3 

^02-^12 ^ 


ill-l pi pl+^ I p*72-i pi R^+1 
-^10-^20 + ® -^21-^01 


+ 


jpl I jpl I jpl ZT^+l 

-^00-^00 +-^11-^11 +-^22-^^ 




^22-^22 


(4.14) 


This Hamiltonian describes interactions which differ from systems we have previously en¬ 
countered, since here the interactions are site dependent. This behavior shows up naturally 


in a non-commutative theory. In [1^|, it was discussed that the y-deformed SYM corre¬ 
sponds to a form of non-commutative deformation of AA = 4 SYM. 
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Figure 1: Spin chain with three sites. The left graph shows the energy spectrum as a function of 
the phase when q = and h = 0. The right graph shows the spectrum as a function of the 
phase 9, when h = e*® and q = 0. 


If all the 7, are equal, the Hamiltonian above will corresponds to our original Hamil¬ 
tonian ( |3.3| ) with q = 0 and h = e*®. The associated R-matrix is 

R(u) = 

- ue-^^ + E[,E{+^] 

- [El^Elt^ + El,E\+^ + El.E^f] 

+ + E\iE\f + El^Eit^] 

+ (1 - u) . (4.15) 

We have checked explicit that ([4.15| ) satisfies the Yang-Baxter equation. This means that 
the theory is integrable! 


In the rest of this section we will discuss the spectrum when the spin-chain Hamilto¬ 
nian (|3.3| ) is either q-deformed or /i-deformed. Figure || shows the spectrum for a three-site 
spin-chain Hamiltonian. The left graph shows how the energy depends on the phase cj), 
with q = and h = 0. The right graph shows instead how the eigenvalues vary as a 
function 6 , when h = e*® and q = Q. Figure shows the same spectra for a four-site 
spin chain. All graphs contain energies which are the eigenvalues of several states. Highly 
degenerate states are generally a sign of integrability because they reflect a large number 
of symmetries in the theory. 

Let us start by explaining the spectra in Figure |^. When h is zero there is only one sinus 
curve while when q is zero there are three sinus curves. The reason is the transformation 
( |4.11 ), since it maps q = and /i = 0 into h = e*® and q = 0 with the relation of the 
phases 9 = — (j) 2Trnj3. Therefore, for each value of q there exist several values of h 

which differ by a phase 2it/3. For (7 = 0, there is a state, independent of the phase, with 
energy three. This state is absent for h = 0. One example of such a state is |000) — |111). 
The “inverse” transformation, see (4.12), of this state is |120) — |201), which is zero due to 
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Figure 2: Spin chain with four sites. The left graph shows the energy spectrum as a function of 
the phase when q = and h = 0. The right graph shows the spectrum as a function of the 
phase 9, when h = e*® and q = 0. 


periodicity. 

The four-site spin chain (see Figure differs substantially from the spin chain with 
three sites. The case q = 0 is completely phase-independent. The reason is the boundary 
conditions. Actually, spin chains with the number of sites differing from multiples of three 
will have spectra which do not depend on the phase. It will just coincide with the spectra 
for the case q = and h = 0. Starting with the case q equal to a root of unity it is 

possible to make a transformation, changing the boundary conditions, such that the phase 
of q is removed |17|. The change in the boundary conditions is then 


|0°)o = |0°)l 


11 o\ _ „i<I> 11 o\ 

|I /o — e \i )l 


and 


|2°)n = e 


— 2$ I oO\ 


(4.16) 


where $ is a phase factor, the exact form of which is not important for our purposes. The 
effect ( 4.16| ) has on the boundary conditions ( 4.13 ) is, when L is one modulo three. 


| 0 " 


= 12 ’^ 


and 


')o = e 


—i<I>| 


)L ■ 


(4.17) 


If we make the shift I V 


—> ggg corresponds to the boundary 

conditions (4.13). The same procedure can be made when L is two modulo three. This 
means that any q equal to root of unity^ can be mapped to any h with the phase 9 = 
TT -|- 2'K'pln -|- 27rm/3. All values of h will then give the same energy spectrum due to 
the fact that p,n and m are arbitrary integer numbers, so the possible values of 9 will in 
principle fill up the real axis. This implies that the energy must be the same for all values 
of 9. For q = g“*27r/3 h = 0 there is a direct map (see (^)) to the case q = 0 and 
h = —g2™/3 which does not change the boundary conditions. The energy spectra for these 
two cases must be the same. Consequently, the spectra for “all” points coincide with the 
spectrum of = g“®27r/3 h = 0. The fact that the shape of the eigenvalue distribution 


’^q = e*'^ is a root of unity iff n(j) = 0mod27r for n an integer. The phase is then (j) = 2'Kpfn where p is 


an integer. 
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changes drastically depending on how many sites there are suggests that a well-defined 
large L-limit does not exist. However, it might still be possible to find a well-defined large 
L-limit if only L-multiples of three is considered. 


5. R-matrix 

We will now try to make a general ansatz for an R-matrix which has the possibility to give 


rise to our Hamiltonian (^). A linear ansatz will turn out to lead to the cases we found 
in the previous section. To find a new solution the ansatz need to be more complicated, 
for instance consisting of hyper-elliptic functions. We are interested in an R-Matrix of the 
following form 

R(u) — dEi^i “t” bEi^i -[- ^ E^^.i 

T cR'iji-t-l ^ -{- cEi-\-\^i ^ 

+ dT'i-|-l,i-|-2 ® T'*,i+2 + dEi^2,i+l <8> T'i+2,i 
+ S.Eij^2,i ® Eij^2,i+l + eR'i,i+2 <8> £'*+1,1+2) 


(5.1) 


where the coefficients are functions of a spectral parameter u 
Written on matrix form the R-matrix is 

/a 0 0 0 

0 


R = 


0 

0 

0 

e 

0 

d 

VO 


b 

0 

c 

0 

0 

0 

0 

e 


0 

b 

0 

d 

0 

c 

0 

0 


c 

0 

b 

0 

0 

0 

0 

d 


0 

0 

d 

0 

a 

0 

e 

0 

0 


e 

0 

0 

0 

0 

b 

0 

c 

0 


0 

0 

c 

0 

e 

0 

b 

0 

0 


d 

0 

0 

0 

0 

c 

0 

b 

0 


0\ 

e 

0 

d 

0 

0 

0 

0 

a/ 


(5.2) 


A natural first step to look for a R-matrix solution is to make a linear ansatz which will 


give the Hamiltonian ( ^) as in ( 3.13 ). 

The Hamiltonian can also be defined through the permuted IH-matrix 

^ = VR, 

where V is the 9x9 permutation matrix. 

If R{u)\u=uo = £ or 9\{u)\u=uo = ^) the Hamiltonian is obtained as 


(5.3) 


H = V-^R{u)\u=uo 


or H = V-^^{u)\u=uo- 


(5.4) 


The linear ansatz below has the property that it gives the Hamiltonian (3.3) in accordance 
with the first formula in 


a{u) = {h*h — k)u + a , 

b{v) = —qu , 

d{u) 

c{u) = {q*q — k)u + a , 

b{u) = —q*u , 

e{u) 

c(u) = (1 — k^u -|- a , 

* 

II 

(i) 

e{u) 
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(5.5) 







with k and a being free parameters, the Yang-Baxter equations turn out to be independent 
of a while they demand k to be k = ^{1 + h*h + q*q). Inserting the linear ansatz in the 
Yang-Baxter equation we find that the equation is satisfied either if 


q = and h = Q or g = 0 and /i = , (5.6) 

where (p and 9 can be any phase, or if the following equations holds 

= (l + , (5.7) 

^-i3<Pr = (l + , (5.8) 

r = ±(l±p), (5.9) 


where we used the notation that q = re*'^ and h = pe*® and let r and p be any real numbers. 
Here we immediately see that the relations between the real parts of q and h are given by 
the last equation, hence we only need to consider which angles are not in contradiction to 
that. The solution is 

q = re*2™/3 ^ h = {I + , (5.10) 

where we once again let r take any real number. Now we would like to see whether there 
exist solutions if an ansatz is made with the permuted version of the R-matrix ansatz (^]^) . 
We obtain 


a{u) = {h*h — k)u + a , 
b{u) = (1 — k)u + a , 
h{u) = {q*q — k)u + a , 


c{p) = -q*u, 
c{u) = —qu , 
d{u) = —qh*u , 


e{u) = hu , 
d{u) = —q*hu , 
e(u) = h*u. 


The conditions from the Yang-Baxter equation read 


q = -q , 


h* = h^ . 


with no restriction on k and a. The only solution to this is 

q ^ _g2^n/3 ^ ^ g2^m/3 ^ 


(5.11) 


(5.12) 


(5.13) 


(or q = Q and /i = 0). This is the other type of solution we expected from the last section. 
The one corresponding to g = — 1 and h = and that one but with twisted boundary 

conditions. Hence a R-matrix with a linear dependence on the spectral parameter u can 
not give us more integrable cases than already found. We need a more general R-matrix 
solution to find new interesting cases. 


5.1 Symmetries revealed 

In order to address the problem of finding the most general solution for the R-matrix ( |5.2| ) 
it is an advantage to make use of the symmetries. We choose the representation 

3 

R = 'Si Si+ oJiSi STi + jiEi S E2i) ■ (5-14) 

i=l 
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All indices in this section are modulo three if not otherwise stated. The generators Si, Ti 
and Ei are 


= 


Ti = 


El = 



So = 


T2 = 


E2 = 



s .= 


T3 = 


E^ = 



.(5.15) 


How the functions in the R-Matrix (^) are expressed in terms of the functions oji, uji and 


7i can be found in Appendix ( A.l ). The generators are related by 


. 2n(l-k) 

SkSi = e-^^^T2k-i 

z27r(i —fc) 

nSi = e - —Ei_k 


SkTi = Ek-i 


■ 27r(l — k) 


EkSi = S, 


k+l 


TkTi = e 
E,T, = 


S2k-l 


l-k 


SkEi = E"-^Sk+i 
TkEi = Tk-i 
EkEi = Ek+i 


(5.16) 


Using these relations it is straightforward to obtain the Yang-Baxter equations which can 
be found in Appendix (|A.2D . A nice feature of these equations is that all of them, except 
the fourth, the fifth and the sixth, can be generated from the first equation through the 
cyclic permutations cvn+i ^n+i 73 and 72 —> > 7i —> >^n +2 ^n+ 2 - The 

remaining three equations are related to each other by the same cyclic permutation. The 
structure of the equations ( |A.2| ) is similar to the Yang-Baxter equations in the eight vertex 
model [Q |3^ 




(5.17) 


for all cyclic permutations (j, k, I, n) of (1,2,3,4). These equations can neatly be represented 
by writing the elements in rectangular objects 


CUn 





£Un 





Wn 

azz 

^k 


azz 

OZn 

CUn 





(5.18) 


Note the beautiful toroidal pattern. The object above should be interpreted as follows. 
The hrst three rows represent the equation (5.17) with the first column representing the 
hrst term in (5.17) 




UnUj'iJ- , 


(5.19) 
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and the next column is equal to the second term in (5.17) 






The next three rows represent another equation of eight vertex model 


(5.20) 



Wn 

^k 




Wn 





£Un 


/ // 


/ // I / // 

UJn^k^^j + ^k^n^l 


= 0 • 


(5.21) 


Our equations can also be represented in terms of similar rectangular objects, with the 
same toroidal pattern 


(5.22) 


U)2 

Oil 

0 J 2 

71 

73 

UI 3 

LVl 

LO 2 

71 

i02 

W 3 

73 

73 

7i 

LO 2 


L 02 

CUi 

71 

73 


UJ2 

oil 

102 

U)2 

0)3 

73 

UJl 

102 

71 


The first three rows give the second equation in ( [A.2| ) with n = 3. The next three rows are 
the seventh equation in (|A.2|) with n = 1. This suggests that the system of equations (A.2) 


should have a nice solution, just like the eight vertex model. The first row determines the 
rest of the entries, thus all equations can be represented with just the upper row. Hence, 
all the 36 equations can be represented by the following rows 


OJn+l 


OOn+1 

71 

73 

(On+2 

(On 

<On+l 

(On 

72 

73 

On+2 

(O 2 

002 

(Ol 

Ol 

03 

O 3 


On+1 


On 

72 

71 

On+1 

O 2 

72 

Ol 

71 

03 

73 

71 

02 

72 

Ol 

73 

03 


The solution to the eight vertex model is a product of theta functions. The cyclicity and 
periodicity properties of the eight vertex model is mirrored into the rectangular object. Due 
to the combination of addition theorems for theta functions and the intrinsic properties 
of the equations, the rectangular objects make it easy to see if an ansatz solves all the 
equations. We believe that the addition theorems for theta functions generating the solution 
of the eight vertex model should be possible to generalize to any even sized rectangular 
object. It would then be interesting to see if those equations are related to an R-matrix of 
arbitrary dimension. 


5.2 A hyperbolic solution 

If the following ansatz oji = , a)j = and 7^ = where we let Qi, Qi and Ki be 

arbitrary constants is made, it leads us to the following solution 


uji = 6 *^ 1 “ , a)i = 6 *^ 2 “ , 7 ^ = ^ 

UJ2 = 6*52“ ^ = gQiw ^ ^2 = e‘3i“ , 

W 3 = ^ (^3 = ^ -y3 = 6*33“ ., 


(5.24) 
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Figure 3: The eigenvalue dependence on the phase for the the Hamiltonian (5.25). 


The following Hamiltonian is obtained from the above R-matrix solution 


H 




+ El 


i 


-^1+1,1+! 




I E^^^ 


E'+'+e- 




'^d+i,j+2 

+ c * d ,.+2 


E'iU + e‘*E‘+2.,H^K+2., 


n/+l 


T7i/ + 1 I ipl 

Ei+i,i+2 + Ei^i 




/+1 


(5.25) 


where + (5ie“*^^/^)/((5f + Qi ~ Q 1 Q 2 ) (we put Q 3 to zero because it does 

not give us any more information). Here we also made use of the fact that the Hamiltonian 
obtained from the procedure (p.l3 ) can be rescaled plus that something proportional to 
the identity matrix can be added. Actually this Hamiltonian can be related with the 


transformation (3.9) to a completely diagonal Hamiltonian, such that it is included in the 
integrable models mentioned in |T^] . In figure |3| the graph to the left shows how the energy 
eigenvalues of the Hamiltonian ( 5.25| ) depends on the phase cf). The graph to the right 
shows the eigenvalues, of the Hamiltonian if we change the sign in front of the second and 
third term in ( 5.25) ), depending on the phase (j). The graph to the right looks very amusing. 
It looks very similar to the graph to the left if that is turned upside down and deformed 
in a considerable symmetrical way. 


6. Broken Z3 x Z3 symmetry 


Relaxing the one-loop finiteness condition (^), by choosing = 0 and qt h 


in the superpotential (2.3) breaks the Z 3 x Z 3 symmetry. The superpotential is 


IFocTr 


o 


( 6 . 1 ) 


where an overall factor is excluded. This superpotential is actually easier to study since 
the dilatation operator has homogeneous vacua |0)|0)... |0) and |2)|2)... |2). The mixing- 
matrix for the anomalous dimensions has the form of a spin-chain Hamiltonian arising 
from R-matrices found by Fateev-Zamolodchikov (or XXZ) [^] and the Izergin-Korepin 
133]. This type of models were considered in |34| even though the authors never completely 
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classified them. They have a [/(l)-symmetry which can be used to get rid of the phase in 
the complex variable h. 

In this setting, there is no longer a cancelation between the fermion loop and the 
scalar self-energy. The additional contribution to the Hamiltonian is of the form (|B.5| ) 
(see Appendix ^ for details). The spin chain obtained from the superpotential is, with 



( 6 . 2 ) 


The term h*h/2, is the fermion loop contribution from the self-energy. We will show that 
for the special values q = —1 and h = e®'^\/2, this Hamiltonian can be obtained from the 
spin-1 XXZ R-matrix. The phase of h is redundant, the energy does not depend on it, 
and can be phased away through the transformation |i) = The R-matrix for the 

XXZ-model is |3| 


R{u) 


/ s 


\ 

t 

r 


T 

a* 

R 

r 

t 


a 

a 

a 


t 

T 

R 

a 

T 


r 

t 

V 


s) 


s = 1 

t = eJ sinh(tt) 
r = J sinh 2r] 

_ ^i(p Tsinhjisinh^ 
sinhfu+r/) 

ry 7 sinh ?7 sinh 277 

sinh('H+? 7 ^ 

rp Tsinh 7 xsinh(bi— 77 ) 

sinh(7i+77) 

a = et + R 
j =_ 1 _ 

sinh(uH-277) 


(6.3) 


where e = ±1. The e in t in (|6.3| ) is added after checking that the R-matrix still satisfies 
Yang-Baxter equation. If we put u = 0, the R-matrix becomes the permutation matrix. 
Thus, a Hamiltonian can be obtained from the R-matrix by the usual procedure H = 
PR'\u=q- Performing the derivatives at the point u = 0 gives 


s' = 0 
R' = - 


t' = e- 


1 


sinh 2 T/ ’ 
cosh rj cosh 2 r/ 


, cosh 277 
r = — 


T' = - 


sinh 2 r/ ’ 
1 


a' = 


1 


sinh 77 ’ 
a' = et' + R! . 


(6.4) 


sinh 77 sinh 2 r] sinh 2 r] 

Multiplying all parameters with sinh 277 , the new variables, evaluated at 77 = vr/d, leads to 

s' = 0, t' = -l, r =0, a=e^^V2, R' = -l, f' = -l, a=0. (6.5) 
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Figure 4: To the left is the spectra for the case = 0 and = h depending on h 

when g = — 1 , and to the right is the spectra for the case all equal to h (up to a constant 
factor) depending on h when g = — 1 


with the corresponding Hamiltonian 


/O 

0 

-1 

±1 

-1 

±1 

0 



0 

e *'^ v /2 


0 

±1 

-1 


-1 


±1 

0 

V 


0 


( 6 . 6 ) 


If we make the choice e = — 1, this is the spin chain Hamiltonian with deformation h = 
e*'^\/2 and g = — 1 ! Looking at the left graph of Figure of a four-site spin chain we see 
that two lines cross at this point. This might, however, just be a coincidence. A special 
feature with q = —1 is that there is a Z 2 -symmetry due to the invariance under exchange 
of the fields <l>o and ^ 2 - 

The right graph shows the same spectrum, but with all couplings equal to h, up 
to a constant factor. The point h = 1 — \/3 is special, since at this point the transformation 
(^) is “self-dual” , which means here that q = q and h = h. 


7. Conclusions 

We have studied the dilation operator, corresponding to the general Leigh-Strassler de¬ 
formation with h non-zero of AA = 4 SYM, in order to find new integrable points in the 
parameter-space of couplings. In particular we have found a relationship between the 7 - 
deformed SYM and a site dependent spin-chain Hamiltonian. When all parameters 7 ^ are 
equal, this relates an entirely g-deformed to an entirely /i-deformed superpotential. For 
g = 0 and the /i = e*®, where 6 is real, we have found a new R-matrix (see EH- 
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We found a way of representing a general ansatz for the R-matrix, with the right form 
to give the dilatation operator, which makes the structure of the Yang-Baxter equations 
clear. The equations can be represented in terms of rectangular objects, which reveals that 
the underlying structure is a generalized version of the structure of the eight-vertex model. 
We presented all values of the parameters q and h for which the spin-chain Hamiltonian can 
be obtained from R-matrices with a linear dependence on the spectral parameter. Most of 
them were related to the g-deformed case through a simple shift of basis with a real phase 
/3, or a shift with a twist with the phase ±27r/3, which reflects the Zs-symmetry. 

We also found a new hyperbolic R-matrix (|5.25 ) which, through a simple change of 
basis, gives a Hamiltonian with only diagonal terms which was included in the cases studied 
in H- We had a brief look at a case with broken Z 3 x Z 3 symmetry and found that the 
matrix of anomalous dimensions can for some special values of the parameters be obtained 
from the Fateev-Zamolodchikov R-matrix. 

We conjecture that the Yang-Baxter equations found for the general R-matrix have a 
solution which is a generalized version of the solution to the eight-vertex model. If this 
solution exists, it is plausible that there will exist more points in the parameter space for 
which the dilatation operator is integrable. To find a general solution to these equations 
would be of interest in its own right. From a mathematical point of view, it is then 
interesting to generalize the solution to an R-matrix of arbitrary dimension. 

The found relationship between the q- and the h-deformed superpotential should be 
visible in the dual string theory, and should also give a clue of what that string theory 
looks like. Another way to approach the problem, as mentioned in Q, is to first find a 
coherent state spin chain and from that reconstruct the dual geometry. The coherent state 
spin chain |^] is valid for small /3, i.e. g ~ 1. We believe that making use of the basis 
transformation (|3.9|) makes it possible to create a coherent state spin chain for g ~ 1 and 


small h acting with the transformation (3.9) on a g close to one gives a new q close to 
one and a new small h. We also hope that due to the relation between vanishing h and 
vanishing q it is possible to write a coherent sigma model for both q and h close to one. It 
would then be very interesting to find the dual geometry, which corresponds to a further 
away deformation of the AA = 4 SYM. 

One other thing of interest is to extend the analysis to other sectors of the theory and 
to higher loop order. In the /3-deformed case it is possible to argue that the integrability 
holds to higher loop order |^], because the dilatation operator is related with a unitary 
transformation to the case of the usual AA = 4 SYM. In the same way can we argue about 
the /i-deformed case, even though we have to consider the induced effects of the spin chain 
periodicity. 
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A. Yang-Baxter equations for the general case 


The functions in the R-Matrix (5.2) are expressed in terms of the functions Wj, w* and 7 ^ as 


a{u) = 7 i(u) + 72 ('u) + 73 (w), 

b{u) = 7i(^/)e'2-/3 + 72(u)e-'2-/3 + ^ 

b{u) = + 7 i(«)e-' 2-/3 + 73 (u) 

c{u) = UJl{u) + 0 J 2 {u) + 0 J 2 ,{u) , 

c(u) = £Di(n) + a) 2 (n) + £D 3 (n), (A.l) 

d{u) = a; 2 (u)e*^""/^ + + u)^{u ), 

d{u) = a)i(u)e*2^/3 + W 2 (u)e-* 2^/3 + a; 3 (u), 

e(u) = + a; 2 (u)e“*^""/^ + a; 3 (n), 

e{u) = a) 2 (u)e*^""/^ + + u) 2 ,{u), ■ 


Yang-Baxter equations from the R-matrix ansatz (5.14) read 


^n+l^n+2^^ ^«.+2^n+l5^2 TS^n^n-l-l T^ 2 ^ri-|-l^n ^n+l'y2^n+l ^n'ys^n+2 

<^n+l^n+2'yi ~ ^n+2^'n+ll'i + 7l^n+2‘^n “ l'i^ndJn+2 + ‘^n73'^n-|-l “ dJn+2lWn+2 
<^n+lw(j+272 “ ‘^n+2i^n+l7r + 72<^n+l<^n+2 “ 7l‘^n+2w"+i + U)n+2lWn+l ~ '^n+ll2^n+2 


^l^'n+1^2 


^l'-^2n+l^'i + ‘^ 2 <^n+ 2‘^0 “ ^2^2n-l'^l + ~ ^0^2n‘^2 


72t^n+l7l + 73 ^^Li 72 + 7l^^n73 “ ‘^l7>2 “ ‘^27n+l‘^3 “ ^^37^1^^! 
7lWn-l72 + 72 W;+i 73 + 73^^71 - “ ^27n+1^3 “ ^37n^l 

^n-\-l^n-\-2^n-\-l ^'n-\-2^n-\-l^n ‘^'n.TsTl “^n+lTlTs Tl^n+l^n+2 T3‘^ri^nH-l 

LDn+2d^n+l^^n+2 “ ^n+li^n+2^n “ ^nl3l2 + ‘^n-H27273 “ 72w(,+2‘^n-Hl + 73‘^nW"+2 
‘^n+ldJn+2^'U2 “ Wn-H2<+lw"+i - UJn+llW-l + ‘^n-H2727l + ll^^'n+l^n+l “ l2^'n+2^n+2 
U:n+2Ul'n(^n+l “ ^nUl'n+2^n + ^n+n'zl'i “ + ll^Wh+2 “ 

^n^^'n+ld^n ~ ^n+l^^'nd^n+2 ~ dJn+2l3l2 + ‘^n7273 “ 72Wn‘^n+l + l'id^'n+2'^'n 
^n+2^ri-\-l^n-\-l ^'f^-\-^^n-\-2^n-\-2 4” ^?^-l-l727l ^?^-l-27l72 4” 7l^ri-l-2^n-l-2 72^n-l-l^n-l-l 


= 0 , 
= 0 , 
= 0 , 
= 0 , 
= 0 , 
= 0 , 
= 0 , 
= 0 , 
= 0 , 
= 0 , 
= 0 , 
= 0 , 
(A. 2 ) 


Here, we have defined ui = a;(u — v), uj' = uj{u) and u)" = u}{y). 


B. Self-energy with broken Z3 x Z3 symmetry 

We will follow the prescription of to compute the contribution to the Hamiltonian from 
the superpotential (|6.1|), when conformal invariance is broken. The additional terms are 
coming from the self-energy fermion loop. 

The scalar self-energy of the vertices is, in Y = 4 SYM, 
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where L = logx (1/e + 7 + logyr + 2). The scalar-vector contribution to this is 
and the fermion loop contribution is • Half of the fermion contribution comes from 

the superpotential; this is the part which will be altered by the extra /i-dependent part of 
the superpotential. Hence, the additional term to the new spin chain, besides the F-term 
scalar part, is 

h*h g‘^{L + l) 


-N : Tr (^p/)i) : . 


(B.2) 


1 q*q Stt^ 

Then, we will have an effective scalar interaction which just comes from the F-term (since 
we have the same cancelation as in the AA = 4 SYM) 


9ym^ . 


(1 -|- q*q) IGvr^ 


: Vp:, 


(B.3) 


where 


+ Tr [qQ*— Qh*— q*h(l)i(l)i(j)2^o] 

+ Tr [h4>i4>iM2 + h*(f) 2 Mi^i + hh*(l)i(t)i^i^i]) . (B.4) 

The plus-minus sign in ( p.3| ) depends on which sign we choose for the superpotential. Since 
all terms are multiplied by the same divergent factor we can set L = logx“^, just as in the 
case of AA = 4. The contribution from the self-energy to the dilatation operator is 

^ ^ {Ell® I + 1® Ell), (B.5) 

1 -I- q q 

and the F-term scalar interaction contribute with 


(d.-eSli+i - 9®.'+wdSi - 


E: 


/ + 1 


V (1 + 9*9) 

+ ) • (B-6) 

We will now consider the case when g = — 1. The total dilatation operator simplifies to 


H = 


/o 

1 h*h 

^ 2 

1 

1 

-qh* 

\ 

1 

1 

4 Wh 

2 


h 

0 

h 


1 h*h 

^ 2 

1 

1 

h* 

1 


1 

1 h*h 

^ 2 

V 


0/ 


(B.7) 


Here we have chosen a relative minus sign between the contribution from the fermion 
loop and the scalar interaction term. 
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